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Abstract 

Multi-dimensional particle-in-cell simulations are used to study the generation of electrostatic 
shocks in plasma and the reflection of background ions to produce high-quality and high-energy 
ion beams. Electrostatic shocks are driven by the interaction of two plasmas with different density 
and/or relative drift velocity The energy and number of ions reflected by the shock increase 
with increasing density ratio and relative drift velocity between the two interacting plasmas. It is 
shown that the interaction of intense lasers with tailored near-critical density plasmas allows for 
the efficient heating of the plasma electrons and steepening of the plasma profile at the critical 
density interface, leading to the generation of high-velocity shock structures and high-energy ion 
beams. Our results indicate that high-quality 200 MeV shock-accelerated ion beams required for 
medical applications may be obtained with current laser systems. 
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I. INTRODUCTION 



CoUisionless shocks are pervasive in space and astrophysical plasmas, from the Earth's 
bow shock to Gamma Ray Bursters, and are known to be efficient particle accelerators [1, 2j, 
even though the details of the acceleration physics are not yet fully understood. The fast 
progress in laser technology is bringing the study of near-relativistic coUisionless shocks into 
the realm of laboratory plasmas. Intense (/ > 10^^ Wcm~^) laser-plasma interactions allow 
for efficient heating and compression of matter [3j and for the generation of relativistic ffows 
relevant to the study of astrophysical coUisionless shocks [4]. 

Apart from the importance of a better understanding of the fundamental physics associ- 
ated with the formation of coUisionless shocks, there has been a growing interest in exploring 
laser-driven shocks as compact particle accelerators [5-9j. Electrostatic shocks can act as 
a "moving wall" as they propagate through the plasma, reffecting background ions to up 
to twice the shock velocity due to the strong electric field associated with the shock front. 
Previous numerical studies of laser-driven electrostatic shocks have shown that the interplay 
between shock acceleration and target normal sheath acceleration (TNSA [ID]), can lead to 
the generation of energetic ions with a broad spectrum [SJ [6l |9] . 

Energetic ion beams from compact laser-produced plasmas have potential applications 
in many fields of science and medicine, such as radiotherapy [111 IE], isotope generation 
for medical applications [13j, proton radiography [H], and fast ignition of fusion targets 
[15]. However, producing focusable, narrow energy spread, energetic beams has proved to 
be challenging. In particular, radiotherapy requires energy spreads of 1-10% FWHM and 
beam energies in the range of 100 — 300 MeV/a.m.u. [16]. 

Recent experimental [17j and numerical [18j results have shown the possibility of using 
tailored near-critical density plasmas to control the sheath fields at the rear side of the plasma 
and generate shock-accelerated, high-quality ion beams. An exponentially decreasing plasma 
profile at the rear side of the target leads to a uniform and low-amplitude sheath electric 
field from the expansion of hot electrons into vacuum [19j. The slowly expanding ions are 
then refiected by the high-velocity shock formed as a result of the laser-plasma interaction, 
leading to the formation of a energetic beam with narrow energy spread [18j. 

In this paper, we expand these recent results [18j by analyzing in detail the different 
plasma conditions that lead to the formation of electrostatic shocks in plasma and their 
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influence in the properties of the reflected ion beams. We consider both the case of ideahzed 
semi-inflnite plasmas with arbitrary density, temperature, and velocity, and the case of laser- 
driven near-critical density laboratory plasmas. We show that electrostatic shocks can be 
formed in strongly heated plasmas by the interaction of two regions of different density 
and/or different velocity, and that ion reflection will occur either for large density ratios or 
for a limited range of relative drift velocities. We then focus on the possibility of driving 
electrostatic shocks in near-critical density plasmas. We show that there is an interplay 
between different physical mechanisms associated with the laser-plasma interaction at near- 
critical density, including laser fllamentation, electron heating, and density steepening. The 
setting up of a fast return current in thin targets is critical to heating the entire plasma 
volume and density steepening plays an important role in launching a shock capable of 
reflecting the slowly expanding background ions. The importance of the plasma scale length 
at the rear side of the target in order to control the quality of the accelerated ion beam is 
also demonstrated. Under optimal conditions, it is shown that this scheme is scalable to the 
production of high-quality (energy spread of ^ 10% FWHM) 100-300 MeV ion beams for 
medical applications with currently available laser systems. 



This paper is organized as follows. In Section [TH we analyze the formation of electrostatic 
shocks and the characteristics of the accelerated ions from the interaction of plasmas with 
different temperatures, densities, and/or flow velocities. We flrst review the theory of shock 
formation and ion acceleration and then use particle-in-cell (PIC) simulations to study the 
properties of the shock and reflected ions as a function of the initial conditions and we 
discuss the possibility of controlling the quality of the accelerated ion beam. In Section 
TTT[ we study the possibility of reaching the required conditions for shock formation and 
high-quality ion acceleration in the laboratory from the interaction of moderately intense 
lasers with tailored near-critical density plasmas. We identify the important mechanisms 
that lead to the formation of a strong shock capable of reflecting background ions and we 
derive the optimal conditions for the generation of high-quality ion beams in laboratory, 
which are validated by multi-dimensional PIC simulations. We explore the scaling of the ion 
energy with laser intensity showing the possibility of generating 200 MeV protons required 



for radiotherapy with current laser systems. Finally, in Section IV, we summarize our results. 



II. ELECTROSTATIC SHOCKS IN PLASMAS 



The interpenetration of coUisionless plasmas of different density, temperature, or velocity, 
leads to a wide range of instabilities and to the formation of nonlinear structures capable of 
trapping and accelerating charged particles. Depending on the exact nature of the instabili- 
ties that mediate these nonlinear structures, different dissipation mechanisms can occur and 
lead to the formation of Shockwaves. Electrostatic shocks are typically associated with the 
excitation of ion acoustic waves (lAW) in plasmas with cold ions and high electron tempera- 
tures. As these waves grow, they start trapping particles, reaching high field amplitudes and 
leading to the formation of a shockwave. If the electrostatic potential energy associated with 
the shock front is higher than the kinetic energy of the upstream ions, these Shockwaves can 
reflect the upstream ions to twice the shock velocity acting as an efficient ion accelerator. 

A. Theory 

To study the formation of electrostatic shocks, we consider the interaction of two adjacent 
plasma slabs with an electron temperature ratio of B = Tgi/Tgo and a density ratio of 
r = Nei/Neo- Electrostatic shock structures can be generated as a result of the expansion 
of plasma 1 (downstream) into plasma (upstream). Here, electrostatic instabilities at the 
edge of the plasmas can develop leading to the build up of the potential at the contact 
discontinuity. Electrostatic shocks can be formed [20l EJJ as dissipation is provided by the 
population of trapped particles behind the shock and, for strong shocks, by the ion reflection 
from the shock front [22j. Kinetic theory can be used to describe such a system, where both 
free and trapped electron populations are taken into account. The ions are treated as a 
fluid. The kinetic theory for the scenario, whereby an electrostatic shock is supported by 
regions/slabs of arbitrary temperature and density ratios has been outlined by Sorasio et. 
al [20] to study the formation of high Mach number shocks. 

The shock transition region is modeled in the reference frame of the shock; the electro- 
static potential increases monotonically from = at x = to = Acf) at X — X ]^ , as 
illustrated in Figure [l] The electron distribution /e(x, Ve) must be a solution of the stationary 
Vlasov equation and can be determined, at a given position x, as a function of the distribu- 
tion of the plasma at the left (xi) and right (xq) boundaries. The free electron population 
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propagating from the upstream to the downstream region is described by a drifting MaxweU- 
Boltzmann (MB) distribution function, with temperature T^o and fluid velocity Vgh (in the 

laboratory frame, the upstream is assumed to be stationary), fefoivo) = ^ ^""^^^ , 

where A^eO is the density of electrons moving from the upstream to the downstream re- 



gion and Vtho is their thermal velocity, deflned as Vtha = ^/^^^ea/^e, with being the 
Boltzmann constant and rrie the electron mass. We assume that the difference between the 
downstream velocity and the shock velocity is much smaller than the electron thermal ve- 
locity and, therefore, that the fluid velocity of the free electrons in the downstream region 
is approximately equal to zero in the shock frame. The free electrons in the downstream 



region have a MB distribution fef i{vi) = - — ^""thi b ei ^ where A^ei is the density 



of electrons moving from the downstream to the upstream region and Vthi is their thermal 
velocity. The trapped electron population is represented by a flat-top distribution function 
feti = ^^""^- ^ following the maximum-density-trapping approximation [21j, which guaran- 



tees fef i{vi = Vc) = feti at the critical velocity Vc = y that discriminates between 
free ('^i < —Vc) and trapped electrons {\vi\ < Vc). The electron velocity at a given point 



follows from energy conservation = y '^o + = ~ y '^i + The electron density 

along the shock transition is calculated by integrating the electron distribution function, 
yielding no{Lp) = A^eoe'^Erfcfy^] for electrons flowing from the upstream to the downstream 
and ni{ip) = iFe^/^Erfcf^ + ~^^eO^\^^/^ for electrons flowing in the opposite 
direction, where ip = ^ ^ and Erfc is the complementary error function. The ion density 
is determined from energy and mass conservation, yielding rii = A^^o/ V^l^^^^^T^^' where 
M = Vsh/cso is the shock Mach number, Cgo = {kBTeo/miy^'^ is the upstream sound speed, 
and rrii and rrie are the ion and electron mass. Using charge neutrality at x = Xq we obtain 
A^eO = A^^o = A^o. 

The ion and electron densities can be combined with Poisson's equation to flnd the 
evolution of the electrostatic potential, which is given by | (^f^) + ^{^) = 0, where x = 
x/Ad, Xd = \J ksTe / 4:7re'^No is the Debye length, and the nonlinear Sagdeev potential [22] 
is given by 

^{if) = Pi{ifi,M) - Pel(^,e,r) - PeoM, (1) 



where Peilv', ©,r) = Pe/i(^,e,r)+Pen(^,e,r) = er/(i+r)(e'^/®Erfcy^+2y^Afe+ 

(8/3)99^/^/ V TT©^ — 1) is the downstream electron pressure, Peo{^,^) = l/(l+r)(e'^Erfcy^+ 



2\J~LpJi\ — 1) is the upstream electron pressure, and Pi{Lp^M) = M^(l — a/Y^^~2lp/AP) is 
the ion pressure. The definition of (f{xo) = and the condition of charge neutrahty at Xq 
impose ^^(99 = 0) = and 1^(99 = 0) = 0, respectively. 

Shock solutions can be found for < 0, allowing for a complete description of the 

shock properties [23j. Ion refiection from the shock front will occur when the electrostatic 
potential across the shock exceeds the kinetic energy of the upstream ions, ecj) > (l/2)m^'^^^, 
which corresponds to the critical value 



Ml 



(2) 



Although ion refiection is not included in this analysis, this critical condition can be used 
to infer the required shock properties, as a function of the plasma parameters, that lead to 
ion refiection from shocks. The critical Mach number, Mc^, for ion refiection can be found 
by solving numerically 



V2Mc^ 



e s^'^Erfc 



Mcr 

V2 



l + ref^ + e^Erfc 



V2Q 



+ 



4M3 



3V27re3 



1 



1 + r 



(3) 



In order to study shock formation and ion acceleration in plasmas where the electron 
temperature is relativistic, we generalize this framework to relativistic temperatures [24j. 
Electrons are described by relativistic Juttner distributions 
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where 70,1 are the relativistic Lorentz factors of upstream and downstream electrons, respec- 
tively, /Xeo = ^eC^/fc^reo, and Ki is the modified Bessel function of the second kind. 

The generalized electron pressures are found by following the same procedure as for the 
non relativistic case and are given by 



Peo(9^,r,/Xeo) 



1 + r 



MeO 



-Me 07 




(7) 



6 



re 
i + r 



^eO 



eiri[//eo/e] A 



^76" 



7+ 



-1-1 



Pen('/',€),r,jUeo) = 



r fj,eoe~ 



(aVa^ - 1 - Log[a + - 1]) 



(8) 
(9) 



i + ri^i[/ieo/e] 

where a = 1 + <^/Heo- In the relativistic hmit, fXeo ^ 1, and we get Pei{(p,Q,T, Heo) = 

^r[^eo(i - ^/e) + + e]/[(i + r)e] and Peo{^,e,r,^,o) = ^{i - ^eo)/(i + r). The 

critical Mach number is given by 
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In the limit of large density ratios (F ^ 1), 



^20 and in the limit of low density 



ratios (F <C 1), Mcr oc ^G/F. Ion reflection can therefore occur for low/moderate Mach 
number shocks provided that F 1 and B ^ 1. 



B. Shock formation 



In order to validate the theoretical predictions for the electrostatic shock structure and 
the conditions for ion reflection we have performed 2D OSIRIS |25j simulations of the inter- 
action of two plasmas with different densities, temperatures, and relative velocity. FuU-PIC 
simulations allow us to understand in a detailed and fully self-consistent way the formation 
of the shock structure and the properties of the reflected ions, as they capture the diflFerent 
kinetic processes involved. 

We model the interaction of two semi-inflnite plasmas and we vary their initial relative 
temperature, density, and/or drift velocity. We consider plasmas with non-relativistic (1 
keV) and relativistic (1.5 MeV) electron temperatures. We use a simulation box with 4098 x 
128{c/(jjpiY ^ where ujpi = \J ^i\n\e^ jm^ is the electron plasma frequency of the denser plasma 
(slab 1), which is located on the left-hand side of the simulation box, between xi = and 
X\ = 204:8c /ujpi. The plasma slab is located between Xi = 204:8c/ ujpi and Xi = 4096c /ujpi^ 
and, therefore, the contact point of the two slabs is at Xi = 2048c/ ujpi. In simulations 
with diflFerent density, temperature, and/or drift velocity between the two slabs, slab 1 is 
always the slab with higher density, temperature, and/or drift velocity, and will correspond 
to the downstream plasma once a shock is formed. Slab thus corresponds to the upstream 
plasma region. The size of the numerical grid is chosen in order to resolve the smallest of 

7 



the relevant plasma scales (either the Debye length or the electron skin depth) with at least 
2 points in each direction. For instance, for Tg = 1.5 MeV, Axi = Ax2 = 0.5c/ ujpi = O.SXd 
and At = 0.3uj~^. We use 9 — 36 particles per cell per species together with cubic particle 
shapes and current smoothing for good accuracy. 

Figure [2] illustrates the ion phase space for different initial density ratios F = 2 — 100 
between the two plasma slabs. For very small density ratios (F 2) the expansion of the 
denser plasma into the more rarefied one drives a nonlinear lAW but no ion refiection is 
present (Fig. [2] a). As F increases, the amplitude of the JAW increases and ion trapping 
becomes evident. Around F = 4, the electrostatic field associated with the leading edge of 
the lAW gets high enough to start refiecting ions from the background plasma (Fig. [2]b). 
For very high F, ion refiection becomes dominant, with the majority of the upstream ions 
being refiected by the shock structure and the trapped component becomes less noticeable 
(Fig. [2]c-e). Both the shock Mach number Msh and the fraction of upstream ions refiected 
by the shock rir/riQ increase with the density ratio F as shown in Figure [3j For a plasma 
with initial constant electron temperature (0 = 1) and no drift velocity, the maximum Mach 
number is observed to be between 1.6 and 1.8. 

We have also studied the infiuence of an initial relative drift between the two plasma 
slabs for F = B = 1 (see Fig. For low relative drift velocities a nonlinear lAW is again 
formed but does not allow for significant particle trapping and no ion refiection is observed 
(Fig. [4] a and b), as in the case of low F. As the relative drift velocity is increased, the 
amplitude of the lAW becomes larger and particles are trapped and refiected by the shock. 
The shock is observed to start refiecting ions for a relative Mach number between the two 
slabs Mi^o = ^ 3 (Fig. [i] c). Again, both the Mach number of the generated 

shock and the fraction of refiected ions increase with the relative drift velocity between the 
two plasma slabs. For F = 1 and = 1, the maximum Msh reached is between 2 — 3 for 
Mi^o — 4 — 5 (Fig. [5]). For very large relative fiows (Mi^o > 10), as the relative drift velocity 
starts approaching the electron thermal velocity, '^1^0 ^ '^t/i, the kinetic energy of the fiow 
is much larger than the electrostatic energy at the contact discontinuity and the fiows are 
only weakly perturbed. For the simulated times {t < lO^cj"^) no shock is formed (Fig. |4]d 
and e). In the opposite limit, when Vi^o ^ Vth^ two-stream and Weibel-type instabilities [4] 
are expected to dominate the shock formation process. 

As the temperature ratio between the two slabs is increased, larger shock Mach numbers 
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can be reached and a wider range of relative drift velocities can lead to the formation of 
electrostatic shocks. For instance, for B = 10 and Mi^o = 10 a shock is formed with Msh = 7 
and for B = 100 and Mi^o ^ 20 shock Mach numbers as high as 20 can be reached. The 
laboratory study of such high Mach numbers [26j would provide important insight on the 
formation of electrostatic shocks in space with Msh = 20 — 100. In simulations where the 
two plasma slabs have different temperatures but the same density and no initial relative 
velocity, no shock is expected and none has been observed. 

C. Ion acceleration 

From the study of the formation of electrostatic shocks for different relative densities, 
temperatures, and drift velocities it is possible to infer the critical Mach number for ion 
reflection, M^r- For a given combination of initial density ratio F and temperature ratio B, 
we have varied the initial drift velocity between the two plasma slabs in order to determine 
the lowest Mach number for which ion reflection is observed, which corresponds to M^r- 
Figure [6] illustrates M^r as a function of F and B. We observe that the critical Mach number 
for ion reflection is in good agreement with theory (Eqs. ^ and ([lO])), as indicated by the 
red and blue circles and crosses in Figure [6| At high density ratios F > 4, the expansion of 
the two slabs (initially at rest) is sufficient to form the shock and reflect the ions. At lower 
density ratios, the plasma slabs need to have an initial relative drift in order to reach M^r for 
ion reflection. The Mach numbers observed in PIC simulations when ion reflection occurs 
lie very near to the theoretical curve for Mcr(F, B) and therefore we can consider that the 
ion velocity will be given by Vions ^ "^McrCgo- The acceleration of ions to high energies in 
the shock requires strong electron heating in order to increase the sound speed. 

In more realistic plasma configurations, where finite plasma slabs are considered, it is 
important to address the expansion of hot electrons into vacuum and the role of the resulting 
space-charge field on the quality of the shock accelerated ion beam. This TNSA field will 
accelerate the upstream ions to a given velocity '^o- The shock will then reflect the upstream 
ions to a velocity Vions — "^McrCgo + '^o- To investigate the role of competing fields in finite 
size plasmas we have preformed 2D simulations where each plasma slab has a thickness of 
200c/ ujpi and are followed by a vacuum region. In the first case, we use a density ratio 
F = 2 (Figs. [T] a and b) and in the second case F = 10 (Figs. [7] c and d). In both cases 
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9 = 1 (Tg = 1.5 MeV). For the abrupt plasma- vacuum transition, the electrostatic field in 
the sheath at the rear side of the upstream plasma is nonuniform and introduces a chirp in 
'Uo, broadening the ion energy spectrum as typical of TNSA [10] (Figs. [7] b and d). This 
sheath field can be controlled by using an expanded plasma profile in the upstream slab. 
For an exponential plasma profile with scale length L^, the sheath electric field is constant 
at early times {t <C 4:Lg/cso) [19j and its amplitude is given by 

Etnsa = (11) 
eLg 

A uniform sheath field will preserve the monoenergetic ion distribution as particles are 
refiected by the shock. This is illustrated in Figs. [7] e and f, where we replace the low 
density slab of Figure [7| a with an exponentially decreasing profile starting from the same 
peak density. The TNSA field is now approximately uniform (Fig. [t] e) leading to a slow 
expansion at uniform velocity of the upstream ions (Fig. [7|f). These expanding ions are 
then refiected by the electrostatic shock and cross the sheath region while preserving their 
narrow energy spread (Fig. [7|f), thus indicating a configuration suitable for the generation 
of monoenergetic ion beams. 

These results indicate that high energy and high quality ion beams can be produced 
from Shockwave acceleration in heated plasmas with an exponentially decreasing density 
profile. In order to achieve good quality in the accelerated ion beam it is necessary to 
guarantee that the velocity of the expanding upstream ions, '^o, is significantly smaller 
than the shock velocity by the time the shock is formed and starts refiecting the upstream 
ions, r^, i.e. Vgh ^ ^lo^r/Lg. For strong shocks, where ion refiection is the dominant 
dissipation mechanism, the ion refiection time, r^, is similar to the shock formation time 
and corresponds to the time an upstream ion takes to accelerate to Vgh in the presence of the 
shock electrostatic field. Viewed another way, in the shock frame, where the upstream ions 
are moving towards the shock at —Vgh^ refiection occurs when the electric field associated 
with the shock is able to stop the incoming ions. For the sake of simplicity, we assume the 
upstream ions initially at rest {vq = 0) and a uniform electric field, Egh = — within 
the shock transition region, L^/^, which for electrostatic shocks is of the order of the Debye 
length, Xd. Let us use Lg^ = ^^A^, with 5 = 0(1). The refiection time can then be estimated 
as Tr = SrriiV shXo / {ecf)) . As we have seen, for shocks driven by the interaction of two plasma 
regions with different densities and low or null initial relative drift velocity, the shock Mach 
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number lies near Mcr and therefore we can use ecj) = {l/2)miv'^^^ yielding 



= (12) 



We note that the obtained expression for the ion reflection/shock formation time is consistent 
with the numerical results obtained by Forslund and Shonk , where the shock formation 
time increases approximately linearly with the Msh before reaching the critical Mach number, 
a for Msh ^ 1-5 the shock formation time is 4:71 /upi. The necessary condition for the 
generation of monoenergetic ion beams can then be written as Lg ^ ^Lsh- 



III. LASER-DRIVEN ELECTROSTATIC SHOCKS 



The conditions required to drive strong electrostatic shocks and generate monoenergetic 
ion beams can be obtained in practice from the interaction of an intense laser pulse with 
plasma. The rear side exponential profile, similar to that shown in Figure [7|e, can be natu- 
rally formed by ionization/pre-heating of the target and consequent expansion, for instance 
due to the laser pre-pulse or an earlier laser pulse of lower intensity. Previous work on 
electrostatic shock formation from laser-plasma interactions focused mainly on laser-solid 
interactions [5l |^ , where electron heating occurs at the vacuum-plasma surface and then 
rely on collisionless plasma processes to heat up the dense background plasma. In this case, 
very high laser intensities are required in order to heat the high density electrons to MeV 
temperatures. Here, we focus on the use of near-critical density plasmas, for which the laser 
can interact with a significant volume of the target and efficiently heat the electrons. 



A. Laser-plasma interaction at near-critical density 

As an intense laser propagates in a plasma with density varying from undercritical to 
critical, ricr^ it will be partially absorbed by heating up the plasma electrons. Depending 
on the laser intensity, polarization, and target density, different absorption and particle 
acceleration mechanisms can occur. For instance, in the underdense region of the target the 
laser can undergo filamentation [28j, self-focusing [28ll29j, and stimulated Raman scattering 
[30] . As it reaches near-critical densities it will steepen the plasma profile locally [31j and 
heat electrons due to a J x B mechanism [32l [33] . 
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Assuming that the laser interacts with the majority of the plasma electrons, the electron 
temperature, aksT^ = eg, can be roughly estimated by equating the plasma electron energy 
density to the absorbed laser energy density, aa^ricLtargetkBTe = V^^iaser^ where a is 3/2 for 
non-relativist ic plasmas and 3 in the relativist ic case, r] is the absorption efficiency, and the 
relativistically corrected critical density aoUc has been used, yielding 

re[MeV] ^ 0.078^ao-^^^^^4^- (13) 

For relativistic laser intensities, ao > 1, and steep density profiles at the laser-plasma inter- 
action region, the temperature of accelerated electrons is expected to be close to pondero- 
motive [3l [33], which leads to a similar dependence (x Qq. Under these conditions and 
for a typical target size Ltarget < 1 nim, laser pulses with picosecond duration can heat the 
plasma electrons to MeV temperatures, leading to high shock velocities and high refiected 
ion energies. 

In order to drive an electrostatic shock, apart from providing the electron heating, it is 
necessary to create a sharp density variation and/or a relative drift velocity between different 
regions of the plasma as seen in Section [llj The plasma push and density steepening due to 
the radiation pressure can provide the required conditions. As the laser is stopped around 
the critical density and steepens the plasma profile, the heated electrons propagate through 
the back side of the target, where they find unperturbed plasma at a similar density, driving 
a return current that pulls the background electrons to the laser region where they are 
accelerated. Therefore, thin targets with peak density around the critical density allow for 
an efficient heating of the entire plasma. 

The initial build up of the return current together with the quick recirculation of the 
heated electrons due to the space-charge fields at the front and at the back of the target, 
will lead to a uniform temperature profile [6l [34] , which is crucial in order to have a uniform 
shock velocity and a uniform ion refiection. Therefore, the target thickness, Ltarget ^ should be 
limited in order to guarantee that electrons can recirculate in the target before ion refiection 
occurs. For a ion refiection time = 4:71 /ujpi (low Mach number shocks [27j), the limit 
on the maximum target thickness is given by I/target < Syrc/c^p^, or equivalently for critical 
density plasmas Ltarget < \o{mi/meY^^ . 



As noted in the Section II, in order to control the strong space-charge fields and maintain 



a narrow energy spread, it is important to have a large scale length at the rear side of the 
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target. Therefore, the optimal target thickness should be close to the maximum thickness 
for uniform heating. For a symmetric target expansion {Ltarget ^ 2L^), the optimal target 
scale length for uniform electron heating and ion reflection is then [I8j 

B. Shock formation and ion acceleration 

In order to explore the physics of laser-plasma interaction at near-critical density and 
to validate the proposed scheme for the generation high-velocity electrostatic shocks and 
high-quality ion beams, we have performed 2D OSIRIS simulations. In this case we use a 
larger simulation box in order to accommodate a vacuum region on the left hand side of the 
target, where the laser interacts with the plasma, and an extended vacuum region on the 
right hand side, where the plasma will expand and ions will be accelerated. The simulation 
box size is 3840 x 240 {c/ujqY and is resolved with 12288 x 768 cells. 

We model the interaction of a Gaussian laser pulse with a duration of 1885cl;^'^ (FWHM), 
inflnite spot size, and a normalized vector potential ao = 2.5 with a plasma with peak density 
of 2.5nc. The pre-formed electron-proton plasma proflle has a linear rise over IOAq and falls 



exponentially with scale length Lg = 20 Aq (chosen according to Eq. (14)). 

Figure [8] illustrates the temporal evolution of the interaction. At early times, it is possible 
to observe the fllamentation of the laser in the underdense plasma and strong electron heating 
(Fig. [sja and i). As the peak laser intensity reaches the critical density region, there is a clear 
steepening of the local density inside the fllaments where the fleld is amplifled. At this point, 
the peak density is increased by a factor of 2 — 4 and followed by the exponential proflle, 
similar to the case of Figure [7] e, which is critical for the shock to be driven. We note that 
the ions also gain a drift velocity at this critical density region due to the space-charge fleld 
caused by the electron acceleration. This drift velocity is measured to be ^ 0.02c (Fig. [s] 
r), which is slightly smaller than the hole-boring velocity [3j Vhb = clq^/ {ncr/'2np){me/mi) = 
0.026 and corresponds to a Mach number of ^ 0.4 for the measured electron temperature, 
which is 2.2 MeV. Both the density jump and the drift velocity will contribute to the shock 
formation and the interplay between these two eflFects can be controlled by tuning the laser 
and plasma parameters. For the proflle used, and taking into account the results obtained 



in Section we expect the density jump to be the dominant effect in our case, and we 
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observe an electrostatic shock being formed as the result of the expansion of the heated and 
tailored plasma profile (Fig. |8]o). 

Although the majority of the laser light is stopped and cannot interact with the electrons 
at the rear side of the target, a return current is set up due to the current imbalance produced 
by the fast electrons in the unperturbed plasma. The cold electrons at the rear side of the 
target are then dragged towards the laser region where they are heated. In Figure [8] j it is 
possible to distinguish between the population of fast electrons that propagate in the rear 
side of the target and the bulk of the background electrons that have negative momentum 
and are being dragged towards the laser due to the electric field that is set up in the plasma 
(Fig. |8]n). This leads to the heating of the entire plasma volume and, together with the 
electron recirculation provides a uniform temperature as can be seen in Figure [8] k for late 
times. The fraction of laser light absorbed into the plasma is measured to be 60% {rj = 0.6). 
As the uniformly heated plasma expands and a shock is formed, it is also possible to observe 
that the filamented density structures caused by the laser interaction are smeared out and 
the shock front becomes relatively uniform. By this time, the laser interaction is finished, 
and the shock moves at a relatively constant velocity, which is measured to be 0.19 c (Fig. |9]) 
and corresponds to Mgh = 1-7 for the measured upstream temperature T^o = 1.6 MeV. The 
measured Mach number is in good agreement with the theoretical Mcr for large F and B ^ 1, 
Mcr ^ 1.5 — 1.8 (Fig. [6]). The shock structure has a strong localized electric field at the 
shock front, with a measured thickness of Lgh ^ 4:Xd = IOc/uoq^ where = ^JkuT^/ ^T^Upe^ 
is the Debye length, which is much smaller than the mean free path for particle collisions 
{Lsh < K^ - c/v,, - 2 X 10^ A,,, A,, - Cs^/vu - 2 x lO^A^,, for T, = 1 MeV, T, = 100 eV, 
and Tie = rii = 10^^ cm~^). Ahead of the shock, the TNSA field is approximately constant 
and in very good agreement with Eq. ( [TT| ) (Fig. [sjb). The density and field structure is 
similar to the case of Figure [7] e and f , where no laser is used and a denser slab expands into 
a more rarefied one with an exponentially decreasing density profile. 

As the shock moves through the upstream expanding plasma it refiects the fraction of 
the upstream ions which have kinetic energy lower than the electrostatic potential energy of 
the shock to a velocity of 0.26 c (Fig. |9]), which is twice the shock velocity in the upstream 
frame plus the plasma expansion velocity Vq^ producing an ion beam with 31 MeV and an 
energy spread of 12% (Fig. |8]t and Fig. [lO]c). The uniform shock velocity obtained under 
optimal conditions is crucial to get a uniform velocity in the refiected ions as we can see in 
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Figure [9| The reflected ion beam contains approximately 10% of the upstream ions, which 
is consistent with the reflected fraction observed in the interaction of two plasmas with 
moderate density ratios (Fig. [3]). The laser to ion beam energy conversion efliciency is 3 
%. We note that while a high reflection efliciency is desirable in order to accelerate a large 
number of ions it can have a deleterious effect for the beam quality, since, as previously 
noted [35j, the strong dissipation of the shock will lead to a decrease of its velocity and a 
chirp in the ion spectrum. Therefore, moderate reflection efliciencies, which are obtained 
for moderate density ratios/drift velocities, are preferable for the generation of high-quality 
beams. 

We have varied the scale length of the rear side of the plasma in order to validate the 
optimal conditions for the generation of high-quality beams. We observe that for shorter 
scale lengths the TNSA flelds become dominant leading to a very broad spectrum. For 
Lg = Lgo/2 the reflected ions have an average energy of 47 MeV, which is similar to the 



case of a sharp plasma- vacuum transition (Fig. 10 a), but the energy spread was increased 



to 36% (Fig. 10 b). For a larger scale length {Lg = 2Lgo)^ where it is harder to uniformly 
heat the entire plasma, the reflected beam has an energy of 17 MeV and an energy spread 
of 30% (Fig. [lO]d). For a very long scale length (uniform proflle) the laser cannot heat the 
entire plasma region and no shock is observed. 

We have tested the impact of the laser spot size in the shock formation process and in 
ion acceleration. Driving a stable shock front and a stable acceleration requires that the 
shock width (which is close to the laser spot size Wq) is large enough such that the plasma, 
expanding transversely at c^, does not leave the shock width region before the acceleration 
occurs. Assuming an isothermal expansion, this condition yields Wq > Lgo/Mgh^ which for 
Msh ^ 2, gives Wq > IOAq. Simulations performed for the same laser and plasma parameters 
but using a super-Gaussian transverse laser proflle with Wq = 16 Aq, led to the generation 
of a stable shock and a reflected ion beam with 28 MeV and an energy spread of 9 %. The 
energy coupling efliciency from the laser to the ion beam was 2%. Assuming cylindrical 
symmetry, the total number of accelerated ions as inferred from the simulation is given by 
Anions ^ 10'^°(Wo[Mni])^/Ao[//m], where Wq is the laser spot size. This number of ions per 
bunch is ideal for most applications. For instance, in radiotherapy ^10^ ions per bunch are 
used in multi-shot treatment and ^ 10^^ ions per bunch in single shot treatment pT[ [T6] . 

The intrinsic ion beam divergence associated with the shock acceleration process can be 
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estimated if we take into account that the velocity of the accelerated ions in the component 
perpendicular to the shock propagation direction is given by the thermal ion velocity of 



the upstream plasma and the parallel component is given by approximately twice the shock 
velocity. The half-angle divergence is then 9 = tan~^ ' 



J_ Ii 

2M V Te 



. For typical moderate 
Mach numbers (M 2) and electron to ion temperature ratios (in our simulations we 
observe Tg/T^ > 10), we expect an half-angle divergence < 4.5°, which is consistent with the 
observed values of 2° in experiments [17j and 4.1° in simulations where a super-Gaussian 
transverse laser profile has been used. For Gaussian transverse laser profiles the shock front 
will have a larger curvature which will increase the overall beam divergence, since away from 
the laser propagation axis the acceleration will occur at an angle. 

We note that in 3D the dynamics associated with the laser-plasma interaction in the 
front of the target (such as self-focusing and filamentation) and with the formation of the 
space-charge field at the rear side of the target will be different than in 2D. The spot size of 



a self-focusing laser in a plasma is given by Vl^ = VI/o^/l ^^/^o^ where Zq = zr/ ^P/Pc — 1 
is the typical distance for self- focusing, zr = ttWq/Xq is the Rayleigh length, P is the laser 
power, and Pc[GW] = ITnc/n^ is the critical power for self-focusing [36j. For typical high- 
power lasers {P > 10 TW) and underdense plasmas {ric/rip ^ 10), P/Pc ^ 1. For a laser 
spot size capable of driving a stable shock (Wq > IOAq) the typical self-focusing distance is 
then zo > ^i^Ao- This means that it is important to keep the characteristic rise length 
of the plasma profile below a few 10 Aq (which is comparable to the optimal scale length of 
the rear side of the target, Lgo) in order to guarantee that self- focusing does not play an 
important role. On the rear side of the target, the TNSA field amplitude will be smaller 
in 3D which should benefit the generation of high-quality shock- accelerated ion beams. 3D 
PIC simulations of this acceleration process are certainly desirable in order to investigate in 
detail the role of 3D eflFects in the acceleration process. 



C. Ion energy scaling 

It is of great interest to study the potential of shockwave acceleration to generate ions 
in the energy range of 100 — 300 MeV/a.m.u. required for medical applications [16j. As 
the electron temperature increases with increasing laser intensity, it should be possible to 
generate larger shock velocities and high energy ion beams. 
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The final ion energy is given by the contribution of both the shock acceleration and the 
uniform expansion of the upstream plasma. In the relativistic case, the final ion velocity is 
Vions = {Vsh + ^o)/(l + ^^/i^o/c^), where = {2Mcso)/{l + M^c^q/c^) is the velocity of the 
refiected ions in the upstream frame and Vq is the upstream velocity at the shock acceleration 
time tacc- Taylor expanding for c^o/c <C 1, the proton energy for optimal conditions is given 

by 



[MeV] ~ 2M,;reo[MeV] + M, 



tacc (2reo[MeV])3/2 



+ 



PiM!. (15) 

mi/me 



To investigate the ion energy scaling, 2D simulations have been performed for increasing 
laser intensities and the same optimal plasma profile. The peak density was increased 
together with the intensity in order to compensate for the relativistic transparency. The 



electron temperature is observed to scale linearly with the laser amplitude (Fig. 11), which 



agrees with Eq. (13) for a laser to electron coupling efficiency rj = 0.51 (consistent with our 
measured laser absorption). 

For the increased laser intensities increased ion energies are observed up to 512 MeV for 



ao = 20 (Fig. 12). The final energy spread varies between 10% and 25%. The ion energy 



scaling with ao is consistent with Eq. (15) for an acceleration time of tacc = 5500cJo ^ (average 



acceleration time in our simulations). At low intensities the acceleration is dominated by the 



shock refiection (first and second term of Eq. (15)), but at higher intensities the contribution 



from the ion expansion (third term of Eq. (|15|)) also becomes important, leading to a 
transition form a scaling with a^^^ to a^. This favorable scaling allows for the generation 
of high quality ^ 200 MeV proton beams required for medical applications with a 100 TW 
class laser system (ao = 10). 

The generation of 100s MeV ion beams using the proposed scheme can be readily tested 
experimentally at different facilities where laser systems capable of delivering 100 TW to 1 
PW power and pulse durations of 0.5 ps — 1 ps are available. The expanded plasma profiles 
required (10s /xm scale and ^ 10^^ cm~^ peak density) can be obtained from the irradiation 
of a /xm scale solid foil by a first low-intensity laser and subsequent target expansion. The 
use of CO2 laser pulses (Ao = 10 /xm) is an alternative possibility [TTj, allowing for the use 
of gas targets where the required plasma profiles, with mm scales and rie ^ 10^^ cm~^, can 
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be naturally obtained from the ionization of the gas by the laser pre-pulse (or by a train of 
pulses). The use of gas targets has the important advantage of allowing for high repetition 
rates in comparison with the conventional solid targets used in ion acceleration experiments. 



IV. CONCLUSIONS 

We have studied the generation of electrostatic shocks in plasma and the use of these 
shocks to accelerate ions to high energy with low energy spreads. Ion reflection can occur 
for electrostatic shocks driven by the interaction of plasma regions with large density ratios 
or moderate relative drift velocities. The energy and number of the reflected ions increases 
with the density ratio or relative drift velocity. For a finite size plasma, it is important 
to control the sheath field at the plasma-vacuum interface, and that can be achieved by 
having an expanded plasma profile with an exponentially decreasing density gradient. In 
this case, TNSA fields will be approximately uniform and of low amplitude, allowing for a 
slow expansion of the ions that are then refiected by the Shockwave as it reaches the rear 
side of the plasma. 

We have shown that the required conditions to drive strong electrostatic shocks in the 
laboratory can be obtained by interacting an intense laser with a near critical density tailored 
plasma. The laser is absorbed near the critical density interface, leading to a local density 
steepening and heating of the plasma electrons. The fast electrons propagate to the rear 
side of the target driving an electric field due to the current imbalance that drags the 
background electrons from the rear side to the laser region. For thin targets, this allows 
for an efficient heating of the plasma volume. As the heated plasma expands, with an 
exponentially decreasing density profile, the electrostatic shock can refiect the background 
ions leading to the generation of a high-energy and high-quality ion beam. The scale length 
of the plasma profile greatly infiuences the quality of the accelerated particles. 

It was demonstrated that by increasing the peak density of the plasma in order to com- 
pensate for relativistic transparency it is possible to scale this acceleration scheme to the 
generation of 100s MeV ion beams with current laser systems (ao ^ 10), which can have an 
important impact for radiotherapy with compact systems. 
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FIG. 1. Steady state electrostatic shock structure as seen from the shock frame. Electrons from 
the upstream region move freely, while electrons from the downstream region can be either free 
or trapped. Ions, which flow from upstream to downstream, are slowed down by the electrostatic 
potential, and reflected back into the upstream for strong shocks. 
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FIG. 2. Ion phase space structure as a function of the initial density ratio F between two plasma 
slabs/regions for 6 = 1 and Tg = 1.5 MeV. Snapshots are taken at t = 2450 uj~^ . At t = there 
is no relative drift between the two slabs. 
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FIG. 3. Shock Mach number (solid lines) and fraction of ions reflected from the upstream (dashed 
lines) as a function of the initial density ratio F between two plasma slabs/regions for O = 1 and 
^1,0 . 
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FIG. 4. Ion phase space structure as a function of the initial relative drift between two plas 
slabs/regions for F = O = 1. Snapshots are taken at t = 2450 ojp^. 
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FIG. 5. Shock Mach number (solid lines) and fraction of ions reflected from the upstream (dashed 
lines) as a function of the initial Mach number of the relative drift between two plasma slabs/regions 

for r = e = 1. 
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FIG. 6. Critical Mach number for ion reflection in electrostatic shocks as a function of the density 
ratio r and temperature ratio O between the two plasma slabs/regions, for Tgo = 1 keV (dashed 
line [20j) and Tgo = 1-5 MeV (solid line Eq. (p!o|)). The symbols indicate the simulation values 
for the non-relativistic (+) and relativistic (o) electron temperatures, which were obtained by 
measuring the speed of the shock structure (density jump or electrostatic field) when ion reflection 
is observed. 
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FIG. 7. Electric field structure and accelerated ion spectrum from the interaction of two finite 
plasma slabs with a density ratio (a,b) F = 2, (c,d) F = 10, and (e,f) F = 10 followed by an 
exponentially decreasing profile. Initially, O = 1 ( Tg = 1.5MeV) and vi^q = 0. Left panels show 
the initial density profile (black) and early time longitudinal electric field (blue), whereas the right 
panels show the ion phase space (orange) and the spectrum of ions ahead of the shock (black line) 



at t = 7700 uj-^. 
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FIG. 8. Temporal evolution of the laser-plasma interaction at near critical densities, from electron 
heating to shock formation, and ion acceleration. Row 1 shows the evolution of the ion density 
profile and row 2 shows a central lineout of the density along the laser propagation axis. Row 3 
illustrates the evolution of the electron phase-space, row 4 the longitudinal electric field, and row 
5 the ion phase-space. 
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Time [1 / coq] 

FIG. 9. Time evolution of the ion density (green) and longitudinal electric field (orange). The 
strong feature between 3 x IO^ojq^ and 4 x IO^cJq"^ is associated with the laser plasma interaction 
and the fields driven by the fast electrons. The solid line follows the shock and the dotted line 
follows the reflected ions. 
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FIG. 10. Ion phase-space and spectrum shock accelerated ions (dashed hne) for upstream plasmas 
with different scale lengths: a) Lg — (sharp plasma-vacuum transition), b) Lg — I/go/2, c) 
Lg — LgQ, and d) Lg — 2LgQ. The initial density profile is indicated by the solid lines and LgQ is 



given by Eq. (14). 
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FIG. 11. Electron distribution for different laser intensities corresponding to ag = 2.5 (green), 5 
(light blue), 10 (red), 15 (orange), and 20 (blue). The distributions are fitted to a 3D relativistic 
Maxwellian of the form 7(7) = a7^e~^/^^ (dashed lines), b) Scaling of the electron temperature 



with the laser amplitude ag. The obtained scaling is consistent with Eq. (13) for a laser-electrons 
energy coupling efficiency 77 = 0.51. 
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FIG. 12. a) Spectrum of shock accelerated ion beams for different laser intensities corresponding 
to ao = 2.5 (green), 5 (light blue), 10 (red), 15 (orange), and 20 (blue), b) Scaling of ion energy 
with the laser amplitude qq. The obtained scaling is consistent with Eq. ([T5|). 
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